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ABSTRACT 



We present a new approach to constrain the spectral energy distribution of the intergalactic UV background observationally by study- 
ing metal absorption systems. 

We study single-component metal line systems exhibiting various well-measured species. Among the observed transitions at least two 
ratios of ionization stages from the same element are required, e.g. C m/C iv and Si m/Si iv. For each system photoionization models 
are constructed varying the spectrum of the ionizing radiation. The spectral energy distribution can then be constrained by comparing 
the models with the observed column density ratios. 

Extensive tests with artificial absorbers show that the spectrum of the ionizing radiation cannot be reconstructed unambiguously, but it 
is possible to constrain the main characteristics of the spectrum. Furthermore, the resulting physical parameters of the absorber, such 
as ionization parameter, metallicity, and relative abundances, may depend strongly on the adopted ionizing spectrum. Even in case of 
well-fitting models the uncertainties can be as high as ~ 0.5 dex for the ionization parameter and up to ~ 1.5 dex for the metallicity. 
Therefore, it is essential to know the hardness of the UV background when estimating the metallicity of the intergalactic medium. 
Applying the procedure to a small sample of 3 observed single-component metal line systems yields a soft ionizing radiation at z > 2 
and a slightly harder spectrum at z < 2. The resulting energy distributions exhibit strong He n hya re-emission features suggesting 
that reprocessing by intergalactic He n is important. Comparing to UV background spectra from the literature indicates that the recent 
model of Madau & Haardt (2009) including sawtooth modulation due to reprocessing by intergalactic He n with delayed helium 
reionization fits the investigated systems very well. 

Key words, quasars: absorption lines - intergalactic medium - cosmology: observations 



1. Introduction 

■ The intergalactic medium (IGM) observed in QSO absorption 
| spectra at redshifts z < 6 is highly ionized. At redshifts z > 1 

OO ■ photoionization is the dominant ionization process. To inter- 
C^l ' pret quasar absorption spectra and to investigate the physics of 

■ m e IGM ionization corrections have to be applied and a thor- 

■ ough understanding of the ionizing radiation field is required. 
^—1 ' Related key questions are the formation of galaxies and the en- 
t— I , richment of the IGM with heavy elements. For studies of the 

■ enrichment history the spectrum and evolution of the UV back- 
. £h \ ground have to be known to estimate the metallicity of inter- 

• galactic absorbers. It is still under de bate how and when met- 
5_i 1 als are distribu ted into the IGM (e.g. iPorciani & Madaul 120051: 
?3 |[Songaila 2006). However, the recently detected lack of Crv ab- 
sorption systems at z ~ 6 suggests that the IGM has not been pre- 
enrich ed to a high level (Be cker et a n i2009t iRvan -Webe r et alj 
2009). Galaxy formation is also influe nced since the cooling of 
the gas is affected by the metallicity ( Wiersma et al.l 120091) as 



well as by the intensity and spectrum of the ionizing background. 
For instance, photoheating suppresses star forma tion in dwarf 
galaxies (e.g.lEfstathiotil d 19921) : iHoeftet all d2006l) . but see also 
lGnedinll20irjh . 

The highly photoionized stage of the intergalactic material 
is maintained by the UV background which is produced by 
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the radiation of quasars and galaxies filtered by absorption and 
re-emission processes in the IGM. Its spectral energy distribu- 
tion (SED) can be estimated from theoretical calculations based 
on the redshift-dependent properties of the sources and the ab- 
sorbers (e.g. the hardness of the QSO spectra, their number den- 
sity, the column density distribution o f intergalactic absorbers, 
etc.). The UV background model from Haardt & Madaul d200ll 
HM01) is usually adopted in absorption line studies. A particu- 
lar application is estimating the metallicity of the IGM, which is 
important to understand the mechanisms and stellar sources of 
metal enrichment. 

Alternative UV background models from other authors 
dFardal et al.lll998t iFauch er-Gigu ereet al.ll2009h and recent up- 
grading (Mad au & Haardt! 120091) reveal the same characteristic 
spectral shape, which is determined by an underlying power 
law slope of the quasar sources with additional contribution of 
galaxies at low energies. Absorption breaks and re-emission fea- 
tures due to the intergalactic Hi and Hen are superimposed 
on the spectrum. The specifications of these characteristics de- 
pend on the details of the handling of the cosmological radia- 
tive transfer as well as on the adopted properties of the quasars 
and galaxies as the sources of radiation. The spectrum in the 
energy range ~ l-10Ryd is in particular important since the 
species easily accessible in the optical have ionization potentials 
in this range (e.g. Cm 3.5 Ryd, Civ 4.74 Ryd, Sim 2.46 Ryd, 
Si iv 3.32 Ryd, Ovi 10. 15 Ryd), which are, therefore, usually 
adopted to trace the metallicity and the relative abundances of 
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the IGM. However, from investigating intergalactic He n absorp- 
tion, which probes the spectral hardness of the ionizing radia- 
tion, it is known that the UV background i s probably fluctuating 
(e.g. iKriss etail l200lt IShull et all l2004t IFechner et all l2006at 
Fechner & Rentiers 120071 : IShull et al.l l2010). These fluctuations 
in spectral hardness can be attributed to variations of the IGM 
density, effects of radiative transp ort, as well as local s ources 
close to the line of sight. Recently, iBolton & Viell d201 lb found 
from numerical simulation, that spatial inhomogeneities of the 
spectral hardness of the UV background have only little im- 
pact on the ionization balances of Si iv and C iv. This is due 
to the longer mean free path of photons with energies close to 
the Hen ionization edge. However, studies of individual inter- 
galactic metal line systems indicate deviations from the standard 
HM01 spec trum, in particular in the important energy range 3- 
4Rvd (e.g . Agafon ova etaDl2005l l2007t IFechner et al.ll2006bt 
iReimers et al.ll2006l) . 

In order to investigate in detail the spectral energy distribu- 
tion of the UV background metal line systems have to be studied 
by means of photoionization modeling. The ionization fractions 
within an absorber will depend on the spectrum and intensity of 
the ionizing radiation as well as on the density of the absorbing 
material. Therefore, the ionizing spectrum can be reconstructed 
to a given level of confidence for appropriate systems. Suitable 
systems have to exhibit various different species, preferably at 
least two of each observed element to avoid degeneracy with the 
elemental abundance pattern. In principle, this can be done us- 
ing absorption features of C n, C m, C iv, Si n, Si m, Si iv, and 
possibly several low-ionization species at redshift z ~ 2. 

In this paper we introduce a new approach to constrain the 
SED of the UV background by extensive photoionization mod- 
eling of metal line systems. The criteria applied to select ap- 
propriate absorption systems are explained in Sect. [2] Our new 
method is presented in Sect. [3] Its potential and limitations are 
illustrated in Sect. |4]by applying it to artificially constructed test 
absorbers. The selected observed absorption systems are inves- 
tigated in Sect.[5]and the results are discussed and compared to 
the most recent model calculations in Sect. [6] We summarize 
our conclusions in Sect. [7] Throughout the paper abundances are 
given in the notation [ X/Y] = log(X/Y) -los jX/Y)^ with solar 
abundances taken from Asplun d et al.l (|2005). 



2. Selection of appropriate systems 

2.1. Selection criteria 

In order to explore the potential of a new approach to constrain 
the spectral shape of the UV background by modeling metal 
line systems it is important to select appropriate absorption sys- 
tems. Suitable systems should be as simple as possible to avoid 
systematics due to simplifying assumptions about the physical 
properties of the absorber. Therefore, we select systems that (1) 
show absorption in one single component only, i.e. without any 
substructure, (2) are unblended and unsaturated in all important 
transitions, and (3) exhibit various species with at least two ra- 
tios of two different ions from the same element. These selection 
criteria are supposed to reduce implicit assumptions about the 
characteristics of the absorbing material like velocity structure, 
multiple gas phases and elemental abundances. 

(1) single-component systems: Systems exhibiting multiple 
component metal line features usually show strong H i absorp- 
tion which has to be distributed to the individual components. In 
principle, the H i can be assigned to the individual components 
assuming the same metallicity throughout the system. However, 



analyses of damped Lya (DLA) systems probed by gravita- 
tionally lensed sight lines showed that the metallicity can vary 
on small scales (e. g. lLopez et al.l 1999; Lopez & Ellison 2003; 
iRichter et alJl2005h . Selecting single-component systems avoids 
any assumptions about the velocity structure and the distribution 
of metals within the absorbing gas 

(2) unblended and unsaturated absorption features: 
Photoionization models are based on the column densities 
estimated from the observed absorption profiles. In particular, 
the ionization parameter U = n Y /rtu, which is the ratio of the 
ionizing photon to hydrogen density, is fixed by a measured 
column density ratio. To keep the uncertainties of the model 
as low as possible it is important to have exact column density 
estimates. Moreover, we will consider the column density 
measurements to evaluate the tested SEDs, and the constraints 
will be more strict if the error bars of the column densities are 
small. Line parameters can be measured best if the absorption 
features are neither saturated nor blended. 

(3) various species with at least two ratios of two different 
ions from the same element: the ionization parameter of an ob- 
served system can be constrained by comparing a measured col- 
umn density ratio to predictions from photoionization calcula- 
tions adopting a radiation field. For a given normalization of the 
ionizing spectrum the ionization parameter is uniquely related 
to the density of the absorber. In principle, for any spec- 
tral energy distribution a solution can be found, even though 
extremely hard/soft spectra will yield unrealistically low/high 
densities and probably unreasonable metallicity values. Thus, 
the derived physical parameters rely on a reasonable choice of 
the radiation field. If two or more column density ratios are 
available, both have to be reproduced at the same ionization pa- 
rameter. This will be impossible for unrealistic ionizing spectra 
providing a tool to evaluate the plausibility of various SEDs. 
In order to be independent of assumptions about the relative 
abundances, which may deviate from solar in th e IGM (e.g. 
lAguirre et al.l2004 120081: rFechner & Richterll2009l) . the adopted 
ratios have to consider two species of the same element. For the 
systems that will be studied in the following the ratios C m/C iv 
and Si m/Si iv or Si n/Si m/Si iv and Al ii/Al m, respectively, are 
available. 

We have searched in published data for metal absorption 
systems satisfying these selection criteria. Three systems at dif- 
ferent redshifts are selected which are suitable for a feasibility 
study. 

2.2. System atz= 1.7529 towards HE 1347-2457 

HE 1347-2457 (z em = 2.61) has been observed as part 
of the ESQ Large Program "Cosmic evolution of the IGM" 
dBergeron et al.ll2004h with UVES at the VLT. The optical spec- 
trum has a resolution R ss 45 000 and a signal-to-noise ratio of 
S IN ~ 35 - 70. F or details of the data reduction we refer to 
lAracil et al.l(l2004l) . 

The system at z — 1.7529 exhibits one single absorption 
component predominately in low-ionization species. Features of 
Cn, Mgn, Aln, Aim, Sin, Si in, Siiv, and also Crv are iden- 
tified. Column densities are estimated by fitting Doppler pro- 
files. The resulting line parameters are listed in the upper part 
of Table Q] and the observed features together with the modeled 
profiles are presented in Fig.Q] 

Unfortunately, Hi is only detected in Lyo- and the corre- 
sponding feature is saturated. Therefore, the column density esti- 
mate may be wrong leading to an incorrect metallicity. However, 
since the absorption is optically thin (i.e. A^hi ^ 10 17 ctrT 2 ), the 
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Table 1. Fitted line parameters. 



QSO 




ion 


z 


\ogN 


fo(km s" 1 ) 


HE 1347 


-2457 


H T 

1 1 1 


1 7^9Qn^^. 


iJ.^+oO ± U.l 1Z 


17.56 


±0.50 






tt 


1 7^98081 
1 . /JZo7(5 1 


i a ns7 j. n ni a 


6.03 


±0.10 






IV 


1 . / jZoV 1 J 


14.0ZU ± U.Ulo 


7.46 


±0.06 






Mg II 




IZ.olO ± U.UUo 


a 

4.JU 


-t- n c\f\ 






Alii 


1.7529035 


11.600 ±0.013 


6.09 


±0.31 






Aim 


1.7529033 


12.229 ± 0.006 


4.64 


±0.12 






Sin 


1.7529005 


12.910 + 0.018 


4.79 


±0.11 






Sim 


1.7529052 


13.553 +0.142 


6.30 


±0.42 






Si IV 


1.7528871 


13.621 +0.007 


6.29 


± 0.05 


HS 1700 + 6416 


1 1 j 

rl I 




i</ii/i_i_ni nn 
JO. 414 ± U.l /U 


26.49 


± 1.55 






Cm 


2.3799345 


13.495 +0.187 


8.64 


± 1.63 






Civ 


2.3799374 


13.059 ± 0.005 


9.69 


±0.14 






Sim 


2.3799212 


11.689 + 0.035 


6.66 


±0.84 






Si iv 


2.3799238 


11.689 + 0.037 


5.74 


±0.99 


HE 0940 


- 1050 


Hi 


2.8265295 


14.612 ± 0.005 


20.13 


±0.22 






Cm a 


2.8265585 


13.232 + 0.013 


6.35 


±0.00 






Civ 


2.8265585 


13.218 + 0.005 


6.35 


±0.13 






Si hi a 


2.8265616 


11.461+0.018 


6.59 


±0.00 






Si iv 


2.8265616 


12.037 ± 0.020 


6.59 


±0.02 



redshifts and Doppler-parameters of C m and Si m have been set to 
the measured values of the C iv and Si iv profiles, respectively 



modeled column density ratios of the metal ions do not depend 
on the exact value of the H i column density. Thus, our main con- 
clusions will not be affected by this uncertainty. However, we 
should keep it in mind for the interpretation of derived metal- 
licity. The features of C iv and Si m as well as Si iv are mildly 
saturated. But since the estimated Doppler-parameters of these 
species agree well, we have confidence in the column density 
measurements. 



2.3. System atz = 2.3799 towards HS 1700+6416 

The metal line systems in the line of sight towards the QSO 



ous authors (e.g. Voael & Reimers 1995; Kohleretal. 


1996 


Petitiean et al. 


1996: iTrippetal. 11997: ISimcoe et al.1 


2006 



Fechner et al. 2006b). We use the co-added spec trum described 
in lFechner et al.ld200 6b) of the data presented in Song ailald 19981) 
and ISimcoe et all (I2002T) . The resolution is R ~ 38 000 with 
S/N ~ 100 at 4000 A. 

The selected system at z = 2.3799 exhibits features of 
Hi, C m, C iv, Si m, and Si iv. The line parameters derived by 
Doppler-profile fitting are summarized in the middle part of 
Table[TJand the observed features are presented in Fig. [2] Again 
the H i absorption is saturated. The column density estimate is 
based on a simultaneous fit of all available transitions of the 
Lyman series (Lya, /3, y). Features of Si in and Si iv are weak 
but clearly detected. The C in profile is slightly saturated and the 
corresponding column density estimate is the most uncertain one 
in this system. 
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Fig. 1. Profiles of the systems at z — 1.7529 towards HE 1347- 
2457. The histogram-like lines represent the observed data. The 
solid lines indicate the profiles from the line fit. 
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Fig. 2. Profiles of the systems at z = 2.3799 towards 
HS 1700+6416. The histogram-like lines represent the observed 
data. The solid lines indicate the profiles from the line fit. 
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Fig. 3. Profiles of the systems at z = 2.8266 towards HE 0940- 
1050. The histogram-like lines represent the observed data. The 
solid lines indicate the profiles from the line fit. 



2.4. System atz = 2.8266 towards HE 0940-1050 

As HE 1347-2457 the QSO HE 0940-1050 (z em = 3.08) is part 
of the ESO Large Program sample. Data quality and reduction 
process are equivalent for both lines of sight. 

We select a system at z — 2.8266 showing absorption in C m, 
Civ, Sim, and Siiv. Again, the features of Sim and Siiv are 
weak and the C in estimate is the most uncertain due to inter- 
lopers from the Lya forest. Redshifts and ^-parameters of C in 
and Si in have been fixed to the values derived for C iv and Si m, 
respectively, to reduce the uncertainties in the column density 
estimates. Using the Lyman series up to Lye the column density 
of H i is well-measured. We summarize the derived line parame- 
ters in the lower part of Tableland present the observed profiles 
in Fig. [3] 

3. Modeling 

Photoionization models are computed for each system using 
CLOUDY fv05.07.06: iFerland et alJ Il99l . The basic idea is 
to construct individual best-fitting models for several spectra of 
the ionizing radiation. Since the systems are selected to provide 
two independent column density ratios, the tested SEDs can then 
be evaluated regarding their consistency with the observations. 
Furthermore, the sensitivity of the resulting physical parameters 
of the absorbers, e.g. density and metallicity, on the spectrum of 
the ionizing radiation can be investigated. 

For a given radiation field the ionization parameter U is de- 
termined by matching the observed column density ratio of two 
species. Because of the uncertainties of the C m measurement we 
adopt Si m/Si iv for all systems. The metallicity of the absorber is 
then adjusted to recover the absolute values of the measured col- 
umn densities. We assume a solar abu ndance pattern and refer- 
ence solar abundances are taken from Aspl und et al.l 12005). The 
tested SEDs are scaled to 7ht = 10 _2L1 erg s _1 cm _2 Hz _1 sr _1 at 
1 Ryd (IScott et al.l l2000). With this normalization the number of 
ionizing photons of each SED is fixed and the derived ionization 
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Fig. 4. Spectral energy distribution of the ionizing radiation in 
comparison to the ionization potential of several species. The 
thick line repr esents the UV backgrou nd at z ~ 2.0 according to 
the model of (lHaardt & Madaull200ll HM01), which is adopted 
as ionizing radiation for test case A (see Sect. 3). The thin line 
indicates an example of a modified spectrum, where the main 
characteristics varied in this work are marked with arrows. They 
are measured in deviation from the HM01 spectrum in dex as 
indicated in the inset. These are (1) the height of the Hen Lya 
re-emission peak at 3 Ryd Aj, (2) the flux level in the energy 
range 3 < E < 4 Ryd A2, and (3) the depth of the He 11 break at 
4 Ryd A3. The presented modified SED can parameterized with 
Ai = +0.5, A 2 = -0.7, and A 3 = -0.5, and is adopted for test 
case B in Sect. [4] 

parameter corresponds to a density «h- Thus, for an individual 
SED a best-fit model is found determined by the resulting values 
of log«H and log U, respectively, the metallicity [Si/H], and the 
relative abundance [Si/C], which is scaled a posteriori to match 
as well the measured column densities of the carbon species as 
accurate as possible. 

This procedure is applied to vario us ionizing spec t ra. At 
first we adopt the UV background from lHaardt & Madaul (|2001 , 
HM01) at the appropriate redshift. Its SED at z ~ 2.0 is indi- 
cated as thick line in Fig. [4] The shape of this spectrum rep- 
resents a mean UV background constituted by the radiation of 
galaxies and quasars which is filtered when propagating through 
the IGM. The main characteristics are a break at 1 Ryd due to 
absorption by intergalactic H 1, a second break at 4 Ryd due to 
absorption by intergalactic He 11, and a small peak at 3 Ryd due 
to He n Lya re-emission. The galaxies included are assumed to 
have an escape fraction of ionizing photons of 10%. Their soft 
radiation mainly contributes to the SED at energi es < 1 Ry d. 

Inspired by the r esults of iFechner et al.l (12006bl) and 
Agaf onova et al.l (120071) . three characteristics of the HM01 spec- 
trum are systematically varied: (1) the height of the 3 Ryd Lya 
peak, (2) the depletion of the intensity in the range 3 < E < 
4 Ryd, and (3) the strength of the 4 Ryd break. We measure each 
of the modifications in departure from the HM01 spectrum in 
dex denoted by Ai, A2, and A3 in the following, as indicated in 
the inset in Fig. [4] The thin line in Fig. [4] shows an example of 
a modified SED with an enhanced 3 Ryd Lya peak by Ai = 0.5, 
intensity depletion at 3-4 Ryd by A2 = -0.7, and a deeper break 
at > 4 Ryd by A 3 = -0.5. 

The parameter space explored is -0.4 < Ai < +1.0 using a 
step size of 0.1 dex, i.e. we vary between spectra with no 3 Ryd 
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case 


logn H 


[Si/H] [Si/C] 


Hi 


Cm 


Civ 


Sim 


Si iv 


CO > 


-3.5 
-3.5 


-1.8 +0.7 
-1.8 +0.7 


15.30 
15.30 


12.659 
13.402 


12.967 
13.215 


11.133 
11.949 


11.512 
12.516 



a absorber ionized by an original HM01 spectrum at z ~ 2 repre- 
sented by the thick line in Fig. [4] 

b absorber ionized by a modified HM01 spectrum with A] = +0.5, 
A2 = -0.7, and A3 = -0.5 represented by the thin line in Fig. [4] 



Lya peak at all (which is at Ai = -0.4 at the probed redshifts) to 
a pronounced peak 1 .0 dex stronger than in the HM01 spectrum. 
The intensity between 3 and 4Ryd is reduced by upto 1.0 dex 
(-1.0 < A 2 < 0.0) with a step size of 0.1 dex. The depth of the 
4Ryd break is modified in the range —1.0 < A3 < +0.4 with a 
step size of 0.2 dex. This means that the intensity at E > 4Ryd 
is either reduced by upto 1 dex or increased by upto 0.4 dex, 
where the latter modification generates a substantially harder 
SED. Since all modifications are combined with each other in 
this feasibility study, 1320 models are computed in total for each 
system. 

The resulting models are evaluated to investigate which of 
the tested SEDs are consistent with the observed metal line sys- 
tems to a given level of confidence. Therefore, we use the second 
observed column density ratio. In case of the systems at z > 2 
this is C m/C iv. Comparing the modeled ratio to the observed 
value log(Ni/Nz) + Alog(A r i/A r 2 ), where Alog(A r 1 /A r 2 ) denotes 
the 1 cr error of the column density ratio, the confidence of the 
best-fit models are estimated. In case of the system at z — 1 .7529 
two more ratios are available (Si 11/Si m and AI 11/ AI in) which are 
applied individually and in a combined manner, respectively, as 
we will present in Sect. [5] Evaluation of the additional column 
density ratio then leads to the SED producing the overall best- 
fitting model. Furthermore, we can derive confidence levels for 
the parameters of the SED modifications. 

4. Test cases 

In order to test the ability of the procedure described above to 
recover the spectrum of the ionizing radiation we construct two 
artificial systems. Presuming an absorber with log«e = -3.5, 
[Si/H] = -1.8, [Si/C] = +0.7, and an observed Hi column 
density of logA^i = 15.3 is exposed to an unmodified HM01 
spectrum at z ~ 2 (corresponding to the thick line in Fig. |4|, 
CLOUDY predicts column densities for C m, C iv, Si m, and Si iv 
listed in Table|2]as case A. As a second test case B the same ab- 
sorber is exposed to a modified HM01 spectrum with Ai = +0.5, 
A 2 = -0.7, and A3 = -0.5, indicated by the thin line in Fig. [4] 
For both test cases we run the same grid of models to establish 
whether the presumed ionizing spectra can be recovered. 

In the following the best-fit photoionization models are 
evaluated considering different aspects about the robustness 
of the method and the influence of the individual parameters 
on the results. To find the best-fit model for every individ- 
ual SED the Sim/Siiv ratio is adopted. The Cm/Civ ratio is 
then used additionally to evaluate the goodness of the model 
with respect to the adopted SED. Assuming an uncertainty of 
AlogAf = 0.01 for each column densities we adopt a 1 cr error 

of Alog(AWjVciv) = V(AlogAfcin) 2 + (AlogAW) 2 = 0.014 
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0.0 1.0 2.0 0.0 1.0 2.0 

Aj - A 2 Aj - A 2 

Fig. 6. Distribution of Ai - A 2 for UV background spectra that fit 
the observed column density ratios up to 1 cr confidence for test 
case A (left panel) and B (right panel). Input values are marked 
with an arrow. 



for the Cm/Civ ratio. Thus, models yielding |(C ni/C iv) Q b s - 
(C in/C iv) mo d| < 0.014 reproduce the column densities within 
1 cr confidence. Here, (C in/C iv) b s denotes the column density 
ratio of the test absorbers given in Table [2] 

4.1. Recovery of the ionizing spectrum 

In Fig. [5]the 1 cr and 3 cr confidence regions of the SED param- 
eters are shown for test case A (upper panels) and B (lower pan- 
els) projected into each of the three possible planes in parameter 
space (Ai versus A 2 , A3 versus A 2 , and Aj versus A3; from left to 
right). Even though the parameters of the input SEDs (marked 
with stars) are within the 1 cr confidence regions, they are not 
recovered unambiguously. 

From the left-most panel it can be seen that the favored Ai 
and A 2 values seem to be correlated. This can be understood 
since both, the height of the 3 Ryd Lya peak and the intensity 
reduction in the range 3-4 Ryd are related to the He 11 recombi- 
nation re-emission. Since the models are based on the Si m/Si iv 
ratio, they are in particular sensitive to the spectral shape at 
3 Ryd since the ionization potentials of Si in (2.46 Ryd) and Si iv 
(3.32 Ryd) are located in this energy range. Thus, the important 
characteristic is the strength of the intensity break at 3 Ryd rather 
than the absolute values of the intensity just below and above this 
energy. Therefore, we introduce the strength of the break a 3 Ryd 
relative to the HM01 spectrum as an alternative parameter, which 
is A] - A 2 . The distribution of Ai - A 2 for the 1 cr-confidence 
models are shown in Fig. [6] For both test cases a clear maximum 
is seen. For test case A we recover (Ai - A 2 ) = 0.2+^4 which 
slightly overestimates the actual value (which is 0.0 since the 
ionizing radiation is the unmodified HM01 spectrum). Test case 
B has (Ai - A 2 ) = 1.2±q*, which is recovered by the procedure 
even though with larger scatter. 

The confidence contours for this parameter versus the 
strength of the intensity break at 4 Ryd A3 are shown in the right 
panels of Fig. [5] The 3 Ryd break A; - A 2 appears to be weakly 
correlated with A3 in the sense that larger A3, i.e. harder spectra, 
favor stronger breaks at 3 Ryd. Since the 4 Ryd break is the He 11 
ionization edge, it is governed as well by the amount of He 11 
in the IGM. Therefore, a correlation between Ai - A 2 and A3 
is expected, since more He 11 means both: stronger He 11 Lyman 
continuum absorption and therefore a stronger 4 Ryd break, as 
well as stronger He 11 recombination re-emission and stronger 
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Fig. 5. 1 cr (marked areas) and 3 cr (dotted lines) confidence regions of the tested SED parameters for test cases A (upper panels) 
and B (lower panels). The three left panels show the projections into each of the three possible planes, i.e. A; versus A2, A3 versus 
A2, and Ai versus A3, respectively. The right panels show the enhancement of the intensity break at 3 Ryd compared to the HM01 
background, which means Ai - A2, versus the strength of the 4 Ryd break A3, where the dashed line indicates a linear fit of the 
correlation (see text). The parameters of the input SEDs are marked with stars. 



Hen L yman series absorption (see Fig. 1 of Mada u~& Haardtl 
(2009) for illustration). Describing this correlation roughly as 
linear we fit Ai - A 2 = (0.50 ± 0.05) • A 3 + (0.22 + 0.02) for case 
A and (0.49 + 0.10) ■ A 3 + (1.40 + 0.05) for case B, indicated as 
dashed lines in the right panels of Fig. [5] In general, the adopted 
method does not lead to reasonable constraints of the intensity 
level at E > 4 Ryd. Models with 1 cr confidence can be found 
for every value of A3 for case A and only the hardest SEDs are 
excluded by the 1 <x constraint for case B. However, a rough es- 
timate of A3 by inserting the recovered A] - A2 into the fitted 
relation yields A3 ^ 0.0 and A3 -0.4 for case A and B, re- 
spectively, which recovers the input value within 0. 1 dex. Fixing 
A3 to the derived value, the best-fit models are then found at 
(Ai,A 2> A 3 ) =* (0.0,0.0,0.0) for case A and (+0.4,-1.0,-0.4) 
for case B, respectively. While the HM01 spectrum is perfectly 
reproduced, the main characteristics of the modified SED (en- 
hanced 3 Ryd Lya peak, depression of the intensity at energies 
> 3 Ryd and > 4 Ryd) are also recovered. Therefore, is should at 
least be possible to give reliable constraints on the spectral shape 
of the UV background also for observed metal line systems. We 
expect, however, to find various SEDs that fit the observed col- 
umn densities to the same level of confidence. Furthermore, we 
should keep in mind, that the parameters of the best-fit models 
are subject to high uncertainties. Strictly speaking, the parame- 
ter ranges marked in Fig. [5] indicate the actual error bars of our 
estimate. 



4.2. Estimated physical parameters 

Usually, photoionization models for metal absorption systems 
are used to derive physical parameters of the absorbing mate- 
rial. Most important are the metallicity and the ionization pa- 
rameter U, which is related to the density «h of the absorber for 
a given normalization of the ionizing spectrum. With our grid 
of best-fit models we are able to investigate the dependence of 
these physical quantities on the shape of the ionizing radiation. 
Since the test cases are based on CLOUDY calculations with 



Table 3. Recovered physical parameters for the test cases. 
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logV 
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0.15 
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-1 33+ 140 
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a recovered SED is the HM01 spectrum, i.e. (Ai,A 2 ,A 3 ) 
(0.0,0.0,0.0) 

b recovered SED has (A,, A 2 , A 3 ) = (+0.4, -1.0, -0.4) 



given ionizing spectra, it is easily verified that the ionization pa- 
rameter (and thus the density), the metallicity, and the [Si/C] 
relative abundance are well reproduced by the CLOUDY mod- 
els if the correct SED is adopted. If system B is modeled with a 
HM01 background as it would be done in a standard analysis, the 
density would be slightly underestimated (by 0.05 dex), but the 
abundances would be heavily overestimated. The HM01 model 
yields [Si/H] = -0.8 and [Si/C] = +1.2 which overestimate the 
actual values by 1.0 and 0.5 dex, respectively. 

Generally, the physical parameters depend on the character- 
istics of the assumed ionizing spectrum. In particular, harder 
SEDs lead to higher metallicities. The marked areas in Fig. [7] 
show the ionization parameter log U (upper panels), the metal- 
licity [Si/H] (middle panels), and the relative abundance [Si/C] 
(lower panels) resulting from the best-fit model of each tested 
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-1.0 -0.5 0.0 

Fig. 7. Recovered physical parameters of the test systems A (left panels) and B (right panels). Marked areas in the upper, middle, 
and lower panels show the ionization parameter log U, metallicity [Si/H] and relative abundance [Si/C], respectively, in dependence 
of the SED parameters Ai, A2, and A3 (from the left to the right) estimated for all tested spectra. The narrower ranges represent the 
results of the models which are consistent with the "observed" column density ratios on a 1 cr level. The actual input values (see 
Table |2| are indicated as stars. 



spectrum versus the SED parameters Ai, A2, and A3 (left, mid- 
dle, right panels), respectively, for test case A (left) and B (right). 
The parameters of the models with 1 cr confidence are indicated 
by the narrower marked areas. Considering all models, harder 
spectra, i.e. larger A3 lead to higher metallicity estimates and 
higher [Si/C] as expected. The ionization parameter depends on 
Ai , i.e. the height of the 3 Ryd Lya peak, in the sense that a more 
pronounced peak leads to lower log U and thus to a higher den- 
sity log «h ■ This is due to using the Si m/Si iv ratio to fix the ion- 
ization parameter since the ionization potential of Sim, i.e. the 
energy to create Si iv from Si in is 2 .46 Ryd, which is affect ed by 
the value of Ai (see Fig. [4] and also Agafon ova et al.l l2007). The 
estimated [Si/C] abundance decreases with increasing 3 Ryd Lya 
peak as well. Virtually no correlation is found between the flux 
level in the range 3 to 4 Ryd and the derived physical parameters. 

Considering all tested spectra the estimated values are spread 
over a wide range (see also Table [3]). Therefore, if the explored 
SEDs represent realistic spectral shapes of the UV background, 
the uncertainties of the derived physical parameters are in the 
order of > 1 dex. The values of the ionization parameter spread 
over A log U = 0.9, where the corresponding densities cover 
Alognn — 1.2. The abundance estimates are scattered over an 
even wider range yielding A[Si/H] = 2.4 and A[Si/C] = 1.8 for 
both test cases. 

The range of possible values can be further constrained if 
only those SEDs are taken into account that fit as well the ad- 
ditional column density ratio C m/C iv within 1 cr confidence 
(dark areas in Fig. [7j. The mean values are listed in Table [3] 
where the quoted error bars refer to the range of covered values. 
The density recovered for case A is slightly underestimated (by 
~ 0.1 dex), while it is slightly overestimated (by ~ 0.2 dex) for 
case B. For both test cases the mean 1 cr-result overestimates the 
metallicity by slightly more than 0. 1 dex. By the same amount 



the [Si/C] abundance is overestimated and underestimated in 
case A and B, respectively. However, within the scatter all values 
are well reproduced. 

The tightest constraints for the physical parameters are ob- 
tained by adopting the best-fit SEDs discussed in the previous 
Section. We have found A3 =2 0.0 and A3 = -0.4 for case A and 
B, respectively. The results obtained when fixing these values 
and selecting the models with 1 cr confidence are also listed in 
Table [3] All parameters are reproduced within 2 cr, however, not 
necessarily to an improved precision compared to the previous 
selection, in particular noticeable for case A. 

For the recovered SED with (A u A 2 , A 3 ) = (0.0, 0.0, 0.0) and 
(+0.4, -1.0, -0.4) for case A and B, respectively, the input val- 
ues are well reproduced. Metallicity and [Si/C] are slightly over- 
estimated (by < 0.2 dex) in case B, which can be explained by 
the slightly harder recovered SED compared to the original one. 

Thus, the shape of the ionizing spectrum as well as the under- 
lying physical quantities like ionization parameter (i.e. density), 
metallicity, and relative abundance [Si/C] of the test cases are 
well recovered by our procedure. We therefore are confident to 
obtain meaningful results also for observed data, provided that 
the SED of the ionizing radiation is, in principle, similar to the 
HM01 spectrum. 



5. Results from observations 

The observed metal line systems presented in Sect. [2] are mod- 
eled in the same way as the test cases where for each system 
the model grid is based on the unmodified HM01 spectrum at 
the appropriate redshift. The ionization parameter of the best- 
fit model for each test SED is determined by matching the ob- 
served Si m/Si iv ratio. In order to estimate the confidence of the 
fit we use in addition C m/C iv for the systems at z = 2.3799 and 
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Fig. 9. Distribution of Ai - A2 for UV background spectra that 
fit the observed column density ratios up to 1 cr confidence for 
the three observed systems at z — 1 .7529 (left panel), z = 2.3799 
(middle panel), and z = 2.8226 (right panel). 



z = 2.8266. For the system at z — 1.7529 both Sin/Sim and 
Aln/Alm are considered (see below). This means that at z > 2 
we probe the same energy range as for the test cases. At z < 2 
the height of the 3 Ryd Lya peak, i.e. A; is expected to have 
the strongest impact on the results since none of the ionization 
thresholds of the species taken into account are at E > 4 Ryd 
(see Fig.|4|l. 



5. 1 . Constraints on the ionizing spectral energy distribution 

The confidence contours with respect to the parameterization of 
the SED projected into all possible planes are presented in Figs. 
[8] and \W\ Obviously, the distributions of the best-fit models of 
the z > 2 systems (Fig. [8} are similar to those found for the 
test cases while the low-redshift system appears to be exposed 
to a different spectrum (Fig. [TOt . This could be due to the dif- 
ferent column density ratios adopted to evaluate the quality of 
the fit but a real evolution of the UV backgroun d is also possi- 
ble dFechner et al.ll2 006b: Ag afonova et al.ll2007l) . Note, that the 
unmodified HM01 spectrum is not included into the confidence 
regions for none of the systems. 



5.1.1. z > 2 

In case of the z > 2 systems the SED parameters of the 1 cr 
models are correlated in a similar way as those of the test cases. 
For both systems more enhanced 3 Ryd Lya peaks, i.e. larger 
Ai are found for less depressed intensity in the 3-4 Ryd range, 
i.e. larger A2. In addition, Ai increases with increasing A3. That 
means harder spectra favor more pronounced 3 Ryd Lya peaks. 

The distribution of Ai - A2, which describes the intensity 
break at 3 Ryd for the most confident models, are presented in 
the middle and right panel of Fig. [9] The mean values are found 
to be <Ai -A 2 > = 0.7±g-| and 0.9±°| for the system at z = 2.3799 
and z = 2.8266, respectively. 

Comparing the intensity break at 3 Ryd Ai - A2 to the break 
at 4 Ryd A3 in the right-most panels of Fig. [8] again yields 
a correlation between those parameters. A linear fit leads to 
Ai - A 2 = (0.44 + 0.04) • A3 + (0.85 + 0.02) for the system 
at z = 2.3799 and Ai - A 2 = (0.31 ± 0.06) ■ A 3 + (1.02 ± 0.03) 
for the system at z = 2.8266, respectively, indicated as dashed 
lines in Fig. [8] Analogously to the test cases, we compute rough 
estimates for the values A3 by inserting the mean Ai - A2 into 
the linear fits. This results into A3 -0.3 for the system at 
z = 2.3799 and A 3 -0.4 for the system at z = 2.8266. 
The most confident models with this 4 Ryd break are found at 



(Ai,A 2 ,A 3 ) = (+0.4,0.0,-0.3) at z = 2.3799, where we have 
interpolated between the models with A3 = -0.2 and -0.4, and 
(+0.3, -0.6, -0.4) at z — 2.8266, respectively. Since the unmodi- 
fied HM01 spectrum has been well reproduced by the test model 
and at least the correct tendencies of the modified UV back- 
ground have been recovered, we are led to the conclusion that 
both systems at z > 2 are ionized by a spectrum different from 
the HM01 background. 

5.1.2. z < 2 

The distribution of the parameters of the most confident models 
for the system at z = 1 .7526 apparently differs from the z > 2 
systems (Fig. [Tol l. Since for this system two additional column 
density ratios are available, we first consider the 1 cr and 3 cr 
confidence regions for the AI 11/ AI m and Si 11/Si m ratios sepa- 
rately (upper and middle panels of Fig. [TO] respectively). Both 
ratios lead to rather tight constraints for the peak at 3 Ryd, Ai, 
since all considered column density ratios (Si ra/Si iv, Si 11/Si m, 
and AI 11/ AI m) are sensitive to the SED at E < 3 Ryd. While 
Aln/Alm requires Ai < +0.4, the Sin/Sim ratio is only re- 
produced with ~ 1 cr confidence by models with Aj = +1.0. 
Furthermore, Sin/Sim favors spectra with A3 > -0.4 while 
AI 11/ AI m leaves A3 unconstrained. 

Since neither the 1 cr nor the 3 cr ranges of AI 11/ AI m and 
Si 11/Si m overlap, we consider a wider range of confident SEDs 
and mark in the lower panels of Fig. [10]those regions in parame- 
ter space that reproduce AI 11/ AI m and Si 11/Si m simultaneously 
to a 5 cr level of confidence. Due to the higher uncertainties of 
the column densities of silicon, the number of confident models 
derived from the Si 11/Si m ratio strongly increases when larger 
cr-ranges are considered. Therefore, the joint 5 cr confidence re- 
gions are mainly determined by the Aln/Alm constraints. The 
3 Ryd Lya peak is then found to be (A] ) = 0.4±Q3. Best-fit val- 
ues of the break at 3 Ryd are in the range (Ai - A 2 ) = 0.9±g g. 
Together with the well-constrained result for Ai we estimate 
A2 ~ -0.5. Furthermore, hard spectra with A3 > -0.2 are 
favored. The actual best-fit model is yielded at (Ai,A 2 , A3) = 
(+0.4, -0.5, +0.4) suggesting that the system at z = 1.7529 fa- 
vors a hard ionizing radiation. 



5.2. Physical parameters 

The estimated values for the ionization parameter log U, the 
metallicity [Si/H], and the relative abundance [Si/C] versus the 
SED parameters are shown in Figs.QT|and[T2]and summarized in 
Table [4] The marked narrow areas in Fig. [TT1 indicate the range 
for those models of the system at z = 1.7529 that fit various 
measured column density ratios simultaneously within 5 cr confi- 
dence, while the narrow areas in Fig.[T2]mark the 1 cr confidence 
regions for the systems at z > 2 based on the C m/C iv ratio. For 
the system at z = 1.7529 (Fig. ITTb these models yield physical 
parameters in an extremely narrow range. On average we find 
(log U) = -2.11 ±0.04 corresponding to (logra H > = -2.34+° °^, 
([Si/H]) = +0.01±°;°j», and ([Si/C]) = -0.86+°;^ where the 
error bars denote the range of estimated values. Adopting the 
HM01 background leads to log U -1.92 and [Si/H] -0.22, 
roughly 0.2 dex higher and lower than the mean of the confident 
models reflecting that those spectra are considerably harder than 
HM01. 

The corresponding plots for the z = 2.3799 and z = 2.8266 
systems are given in the left and right panels of Fig. [12] respec- 
tively. Due to the higher uncertainty of the C in column density a 
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Fig. 8. Projected 1 cr (marked areas) and 3cr (dotted lines) confidence contours of the SED parameters for the observed systems 
at z = 2.3799 (upper panels), and z = 2.8266 (lower panels). The projections in each of the three possible planes are shown (Ai 
versus A2, A3 versus A%, and Aj versus A3; from left to right). The right most panels display the confidence regions considering the 
enhancement of the break at 3 Ryd, i.e. Ai - A2 versus the strength of the 4Ryd break A3, where the dashed lines indicate a fitted 
linear correlation (see text). The contours of the z = 2.8266 system are smoothed for a clearer presentation. 




-1.0 -0.5 0.0 -1.0 -0.5 0.0 -1.0 -0.5 0.0 -1.0 -0.5 0.0 

Ag Ag A 3 A 3 

Fig. 10. Projected confidence contours of the SED parameters for the observed system at z — 1 .7529 considering the AI 11/ AI m ratio 
(upper panels), the Si 11/Si 111 ratio (middle panel), and the combination of both column density ratios (lower panel). The projections 
in each of the three possible planes are shown (A; versus A2, A3 versus A2, and Aj versus A3; from left to right). The outer right 
panels display the confidence contours considering the enhancement of the break at 3 Ryd, i.e. Ai - A2 versus the strength of the 
4 Ryd break A3. The upper and middle panels show the 1 cr (marked areas) and 3<x (dotted lines) ranges of confidence. Since 
there are no SEDs that reproduce the AI 11/ AI m and Si 11/Si m within < 3 cr simultaneously, the lower panels indicate the regions in 
parameter space where both column density ratios are matched within 5 cr confidence 
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Fig. 12. Dependence of physical parameters on the adopted ionizing spectrum for the observed systems at z > 2. The marked 
areas in the upper, middle, and lower panels show the ionization parameter log U, metallicity [Si/H], and relative abundance [Si/C], 
respectively, versus the parameters of the ionizing SED (Ai, A2, A3; from the left to the right) for the best-fit models of the systems 
at z — 2.3799 (left) and z = 2.8266 (right). The narrower areas indicate the parameters of the models consistent with the observed 
C ni/C iv column density ratios on a 1 cr level. 




Table 4. Results for the three studied metal line systems. 



Fig. 11. Dependence of physical parameters on the adopted ion- 
izing spectrum for the observed low-redshift system at z — 
1.7529. The wide and narrow marked areas in the upper, middle, 
and lower panels show the ionization parameter log U, metallic- 
ity [Si/H], and relative abundance [Si/C], respectively, versus the 
parameters of the ionizing SED (Ai, A2, A3; from the left to the 
right) for all best-fit models and those consistent with the ob- 
served column density ratios of AI 11/ AI in and Si n/Si m on a 5 cr 
level, respectively. 
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a recovered SED has (A u A 2 , A 3 ) = (+0.4, -0.5, +0.4) 
b recovered SED has (Ai, A 2 , A 3 ) = (+0.4, 0.0, -0.3) 



c recovered SED has (A,,A 2 ,A 3 ) = (+0.3,-0.6,-0.4) 



wider range of consistent models and parameter values is found. 
For the z = 2.3799 system we derive (log U) = -1.52±9-f2 
for the ionization parameter and a metallicity of ([Si/H]) = 
-2.41±ggj with 1 cr confidence. The results derived for a HM01 
spectrum (log U - -1.50 and [Si/H] -2.12) are within this 
range although the HM01 model is not within the 1 cr confidence 
range. The system at z = 2.8266 yields (log U) = -1.04 ± 0.33 
and ([Si/H]) = —1.31 ±" 70 also cons i stent witn tne HM01 result 
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of log 1/ - -1.00 and [Si/H] -0.73 even though the HM01 
spectrum significantly underestimates the observed C m/C iv ra- 
tio. Note that already the 1 <x results for the metallicity and hy- 
drogen density of the z > 2 systems are spread over more than 
one order of magnitude strongly depending on the hardness of 
the ionizing radiation and the strength of the Lya peak at 3 Ryd. 

Following the analysis procedure applied to the test cases, 
we fix A3 to the recovered value, which is -0.3 and -0.4, respec- 
tively, and select the models fitting at least with 1 cr confidence. 
The ranges of the physical parameters then strongly reduce (see 
Table|4|. In particular, the metallicity is now tightly constrained 
with a spread of A[Si/H] ~ 0.1. The results for the recovered 
SEDs are also listed in Table [4] Compared to the HM01 re- 
sult, we find only slightly different ionization parameters. But 
for both z > 2 systems the metallicity derived from the recov- 
ered SED is lower by ~ 0.3 to 0.7 dex due to the softer spectrum. 
Furthermore, the enhancement of the 3 Ryd Lya peak leads to a 
higher [Si/C] abundance by roughly 0.5 dex. 

6. Discussion 

As shown in the previous Section, evaluating best-fit photoion- 
ization models of three observed metal line systems suggest 
that the ionizing radia tion differs from the usually adopted 
lHaardt & Madaul OOP 1) background spectrum. In two systems 
at z > 2 a softer spectrum has been recovered where the He n 
re-emission peak at 3 Ryd is slightly enhanced and the intensity 
level in the range 3-4 Ryd is reduced. Measuring the deviation 
from the HM01 spectrum in dex, our fiducial SEDs can be pa- 
rameterized with (Ai, A 2 , A 3 ) = (+0.4, 0.0, -0.3) at z = 2.3799 
and (+0.3, -0.6, -0.4) at z — 2.8266, respectively. For one sys- 
tem at z = 1.7529 a hard spectrum even more different from 
HM01 is suggested yielding (+0.4, -0.5, +0.4) as fiducial SED. 

6.1. The spectrum of the UV background 

Our fiducial spectra are shown in Fig. [13] as thick solid lines. 
Keep in mind, that the uncertainties of these SEDs are high and 
hardly to quantify since the parameters are mutually correlated. 
Nevertheless, in the following we compare our best-fit recov- 
ered SEDs to various U V background spectra fro m the literature, 
which are the classical lHaardt & Madaul (120 01) spectrum, vari- 
ous upgradings based on iMadau & Haardt) (2009). a nd the in- 
depen dently computed background fromlFaucher-Gig uere et all 
(2009, see Fig. Q3). 

The recent model of lFaucher- Giguere et al. (2009, FG09) at 
the appropriate redshift is indicated as dotted line in each panel 
of Fig. [13] At each redshift the FG09 background is signifi- 
cantly flatter than our fiducial SEDs and even slightly flatter than 
the pure HM01 background (thin solid black lines). The FG09 
spectra show only very weak recombination emission at 3 Ryd. 
According to the authors this is due to the combined effects of 
the saturation of re-emission in optically thick systems, leakage 
of the re-emitted photons at the ionizing edge, and the frequency 
dependence of the photoionization cross section. In addition to 
the weak 3 Ryd Ly<? peak the FG09 spectra show a high intensity 
in the 3-4 Ryd range, leading to a weak break at 3 Ryd, which is 
even less pronounced than in the HM01 spectrum (thin black 
line in Fig. [T3l ). Since this feature is robustly constrained by 
our models to be larger than in the HM01 case, we conclude 
that the Hen Lyman continuum absorption and recombination 
re-emission may be underestimated in the FG09 models. 

Furthermore, the break at 4 Ryd is weaker in the FG09 
spectra than in the HM01 background. Consequently, the FG09 



spectra provide a har d er rad iation field. As demonstrated by 
IFaucher-Giguere et al] (120091) the depth of the 4 Ryd break 
mainly depends on the slope /3 of the Hi column density dis- 
tribution and the spectral index o-qso of the contributing QSOs. 
The steeper the column density distribution or the softer the QSO 
radiation, the stronger the break at the He 11 edge. In both cases 
the number of He n Lyman limit systems is increased leading 
to more absorption just above 4 Ryd. Thus, comparison to our 
fiducial spectra at z > 2 suggest that either the H 1 column den- 
sity d istributi on slope is stee per than the adopted (3 = 1 .4 of 
Faucher-Gigu ere et al.l (H009) or the contribution of QSOs to the 
UV background at these redshifts is less than assumed. 

The recent finding of a flat, i.e. nearly constant, hydrogen 
photo ionization rate at z > 2 by IFaucher-Giguere et al.l (12008) 
andlDall' Aglio eTail (l2009h suggests that stars contribute signif- 
icantly to the ionizing radiation at z > 3. If galaxies dominate the 
UV background, a softer SED is expected. The studied systems 
may probe the transition from a QSO- to a galaxy-dominated 
background in the redshift range 2 < z < 3. Alternatively, the 
absorbers may be affected by an additional contr i bution from lo - 
cal soft radiation sources (Miralda-Escude 2005; Schave 2006). 

Recently, Mad au & Haardt] J2009) have shown that resonant 
line absorption in the Lyman series of Hen leads to a saw- 
tooth modulation of the UV background spectrum between 3 
and 4 Ryd. The resulting energy distributions at the appropri- 
ate redshifts are also shown in Fig. Q~3] as dashed red lines 
(MH09/HM+S). In this model only QSOs are considered as 
sources of radiation and photoionization equilibrium with a uni- 
form radiation field is assumed. Since the soft radiation con- 
tributed by galaxies is neglected, the sawtooth SED is much 
harder than the HM01 background which takes into account 
galaxy spectra as well. 

A very recent upgrade of the HM spectrum 
(lHaardt & Madaul I20T1I HM11) includes the sawtooth modu- 
lation as well as the radiation of galaxies (for more details we 
refer to lHaardt & Madaul201 ll) . A rough comparison of the new 
model to the HM01 background at z ~ 3 shows that besides the 
sawtooth modulation the HM1 1 spectrum declines slightly more 
steeply between 1 and 3 Ryd and has a weaker break at 4 Ryd. 
At higher energies, the intensity of the HM11 spectrum drops 
significantly, in particular at E > 45 Ryd, while the HMN01 
background is still flat in this energy regime. 

The colored thin soli d lines in Fig. [13] repr esent a saw- 
tooth model according to IMadau & Haardt (2009) which addi- 
tionally assumes late He 11 reionization by artificially boosting 
the Hen/Hi ratio by a factor of 10, 30, and 500 (MH09/DR10, 
MH09/DR30, MH09/DR500), respectively. In these cases, the 
spectrum is suppressed above 4 Ryd and the sawtooth modula- 
tion results in a net depression of the intensity level in the range 
3-4 Ryd. Both effects increase with the amount of existing Hen, 
i.e. with the He n/H 1 ratio. Furthermore, the re-emission peak at 
3 Ryd is enhanced. Thus, at z > 2 the MH09/DR spectra show 
similar overall features as our recovered SEDs. In particular, for 
the system at z = 2.3799 the recovered spectrum near 4 Ryd fol- 
lows closely the MH09/DR30 spectrum (middle panel of Fig. 
[T3l and for the system at z = 2.8266 the 3 Ryd Lya peak of the 
recovered SED and all the MH09/DR spectra are nearly identical 
(lower panel of Fig. IPTb . 

For each UV background spectrum from the literature we 
compute with CLOUDY a best-fit photoionization model for 
each observed system. The resulting physical parameters are 
summarized in Table [5] In addition to the hydrogen density 
log«H, the ionization parameter log U, the metallicity [Si/H], 
and the relative abundance [Si/C], the Hen/Hi ratio log 77 is 
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Table 5. Comparison of UV background spectra from the litera- 
ture. 



-sys 


ionizing spectrum 


log« H 


logt/ 


[Si/H] 


[Si/C] 


log;/ 


a 


1.7529 


recovered SED a 


-2.30 


-2.15 


+0.08 


-0.81 


1.61 


4.7 




HM01 


-2.61 


-1.92 


-0.22 


-0.99 


1.95 


6.8 




FG09 


-2.27 


-2.17 


+0.15 


-0.76 


1.50 


8.5 




MH09/HM + S 


-2.16 


-2.27 


+0.31 


-0.61 


1.29 


12.4 




MH09/DR10 


-2.37 


-2.08 


+0.01 


-0.84 


1.78 


4.2 




MH09/DR30 


-2.44 


-2.01 


-0.09 


-0.91 


1.96 


7.9 




MH09/DR500 


-2.77 


-1.68 


-0.46 


-1.15 


2.51 


15.5 


2.3799 


recovered SED 


-2.77 


-1.67 


-2.39 


+0.46 


2.37 


0.1 




HM01 


-3.02 


-1.50 


-2.12 


-0.11 


2.20 


6.3 




FG09 


-2.71 


-1.72 


-1.81 


-0.17 


1.81 


7.5 




MH09/HM + S 


-2.56 


-1.87 


-1.53 


-0.24 


1.50 


7.9 




MH09/DR10 


-2.88 


-1.57 


-1.81 


-0.22 


2.20 


2.5 




MH09/DR30 


-2.94 


-1.51 


-1.96 


-0.10 


2.47 


0.1 




MH09/DR500 


-2.96 


-1.49 


-2.51 


+0.41 


2.98 


4.3 


2.8266 


recovered SED C 


-3.54 


-0.93 


-1.46 


-0.34 


2.79 


0.1 




HM01 


-3.52 


-1.00 


-0.73 


-0.79 


2.34 


30.9 




FG09 


-3.17 


-1.26 


-0.50 


+0.83 


1.89 


34.7 




MH09/HM + S 


-2.94 


-1.49 


-0.09 


-0.85 


1.53 


31.5 




MH09/DR10 


-3.32 


-1.13 


-0.24 


-1.05 


2.27 


14.5 




MH09/DR30 


-3.43 


-1.02 


-0.33 


-1.05 


2.52 


9.6 




MH09/DR500 


-3.57 


-0.87 


-1.02 


-0.72 


3.13 


3.6 



recovered SED has (A,, A 2 , A 3 ) = (+0.4, -0.5, +0.4) 
recovered SED has (A,, A 2 , A 3 ) = (+0.4, 0.0, -0.3) 

recovere d SED has (A t , A 2 , A 3 ) = (+0.3, -0. 6, -0.4) 

HM01 : lHaardt & Madaul d2001h, FG09 : iFaucher-Giguere et alj 
(2009), MH09/HM+S: iMadau & Haardtl d2009h photoionization 
equilibrium model, M H09/DR10, MH09/DR30, MH09/DR500: 
Madau & Haardt ( 2009) delayed helium reionization model where 
the He n/H i ratio is boosted by a factor 10, 30, and 500, respectively 



listed. The He n/H i ratio measures the hardness of the ioniz- 
ing radiation which is mainly determined by the intensity level 
above 4Ryd. Therefore, the strongly absorbed MH09/DR500 
background yields the highest 77-values indicating a soft radia- 
tion field while the relatively flat MH09/HM+S and FG09 spec- 
tra lead to lower values corresponding to hard spectra. In or- 
der to evaluate the quality of the models we compare the pre- 
dicted column density ratios to the observed values in Fig. [14] 
As above the C 111/C iv ratio is adopted at z > 2 (lower panels) 
while Si 11/Si m and Al 11/ Al m are used at z < 2 (upper panel). 

In case of the low-redshift system the quality of the fit is 
dominated by the Al 11/ Al m ratio (indicated as triangles) which is 
tightly constrained by the observations (narrow marked range). 
Besides our fiducial model, the MH09/DR10 spectrum fits the 
observed column density ratios with 4.2 cr confidence. Though 
our recovered SED suggests a hard ionizing radiation, the mod- 
els based on the harder backgrounds FG09 and MH09/HM+S 
are less accurate. However, the selected 5 cr-confidence mod- 
els cover the range of hardness log 77 = 1 .70±9j 9 correspond- 
ing to 77 « 50±l° . The FG09 and MH09/HM+S spectra are even 
harder yielding 77 ~ 30 and 20, respectively, while the 77-value 
of the MH09/DR10 background (~ 60) matches the range of 
the recovered spectra. Comparing to direct measurements of the 
He n/H 1 ratio at z k 2 (IZheng et alJl2004t IFechner et al.l l2006a: 
iFechner & Reimersl 120071: IShull et all 1201 Oh shows that the re- 
covered 77-values are at the hard end of the expected range. 



A hardening of the U V background at z < 2 has been claimed 
alread y by other studies (IFechner et alJl200"6bl : lAgafonova et al.l 
|2007|) . However, for the system at z = 1.7529 the results should 
be interpreted very cautiously. Absorption features of various 
species are detected in this system (see Fig. Q]). In particular, 
the system exhibits strong lines of C iv and C 11 simultaneously. 
Due to their different ionization potentials (1.79 Ryd for C 11 and 
4.74 Ryd for C iv) both species probably do not live in the same 
gas phase. While the C 11 absorption is thought to arise from the 
same phase as Si 11, Si in, Al 11, and Al m, C rv is expected to trace 
a more highly-ionized gas phase. Possibly, part of the Siiv ab- 
sorption may arise from the C iv phase as well. If so, the mea- 
sured Si m/Si iv ratio used to constrain the models, would actu- 
ally be a lower limit and, therefore, our results would be biased. 
Thus, an independent confirmation for a z < 2 metal line system 
using C m/C iv is desirable. But due to the short restwavelength 
of Cm (977.0201 A) this requires high-quality near-UV obser- 
vations. 

At z > 2 our fiducial SEDs optimally reproduce the observed 
Cm/Civ ratio since they are constructed to do so. Furthermore, 
the delayed reionization sawtooth models lead to surprisingly 
good fits for both systems. The column densities of the system at 
z = 2.3799 are optimally modeled by the MH09/DR30 spectrum. 
The other delayed reionization models lead to acceptable fits as 
well reproducing the C 111/C iv ratio within 2.5 cr (MH09/DR10) 
and 4.7 cr (MH09/DR500), respectively. Note, that the physical 
parameters derived when adopting the MH09/DR30 spectrum 
differ significantly from those estimated with our recovered SED 
by upto 0.5 dex (see Table|5]l. However, the values are consistent 
with the 1 cr range listed in Table [4] The system at z = 2.8266 
is fitted best by the MH09/DR500 spectrum (3.6 <x confidence). 
Though the amount of He 11 required to generate such a soft radi- 
ation field seems unrealistically high, we conclude that the ion- 
izing radiation should indeed be rather soft. 

At leas t for these two indiv idual systems our findings are 
contrary to IVasiliev et alj d2010t) . They have compared the ob- 
served Chi/Civ ratio of 10 absorbers to (non-)equilibrium pho- 
toionization models for several UV background spectra from 
Madau &Haardj {2009). The authors come to the conclusion, 
that the delayed reionization models (MH09/DR) and thus a 
large fraction of He 11 is disfavoured by their data. 



6.2. Physical parameters 

The physical parameters of the analyzed metal line systems are 
also given in Table|5] Though the various adopted ionizing spec- 
tra lead to different results for ionization parameter, i.e. hydro- 
gen density, metallicity [Si/H], and [Si/C] abundance, the de- 
rived values indicate distinct results for the three systems. 

The system at z — 1 .7529 shows roughly solar [Si/H] with 
silicon depleted with respect to carbon with [Si/C] ~ -0.8 for 
all adopted spectra. This means, the system is apparently metal- 
rich with [C/H] ~ +0.8. The system at z = 2.3799 has sub-solar 
metallicity, [Si/H] ~ [C/H] ~ -2.0, with roughly solar [Si/C]. 
The soft recovered SED and the even softer MH09/DR500 spec- 
trum lead to enhanced silicon abundance of TSi/C I ~ +0.4 as 
would be expected in the IGM dAguirre et alJl2008h . The system 
at z = 2.8266 appears to be metal-rich as well. For the adopted 
spectra we find about solar metallicity, [C/H] ~ -0.1, and sil- 
icon being depleted with respect to carbon by [Si/C] ~ -0.6. 
Note, that the hard UV background of FG09 suggests enhanced 
silicon for this system (see Table|5]l. 
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Two of the analyzed sy stems have been investigated in detail 
bv lAgafonova et al.i (120071 z = 2.3799 towards HS 1700+6416) 
and lLevshakov et al.1 d2009l z = 1.7925 towards HE 1347-2457), 
respectively. For both systems the measured column densities 
are g enerally in good a greem ent. In case of the z = 2.3799 sys- 
tem, lAgafonova et alj (120071) also estimate similar abundances 
([C/H] = -1.94 and [Si/C] = +0, 33) even though they applied 
a slightly different UV background spectrum, 

Ho wever, for the system at z — 1 .7529 iLevshakov et al.l 
(2009) estimate a significantly higher Hi column density 
(logAf(Hi) 16.0 in comparison to our ~ 15.5). Since the 
Lya feature is saturated and higher order Lyman series lines are 
not observed, the H i column density cannot be estimated unam- 
biguously. But the difference may explain the different resulting 
metallicitv. ILevshakov etail (120091) derive [C/H] = 0.1 . . .0.2 
and [ Si/C] = -0.3, indicatin g a moderate metal-rich absorber. 
While Levshakov et al. (2009) also apply a hard but slightly dif- 
ferent radiation field, our very high metallicity of [C/H] — h0.8 
may be due to the underestimated H i column density. 

Finding silicon depleted in intergalactic absorption systems 
like for the systems at z = 1 .7529 and z = 2.8266 is rather un- 
usual. According to lAguirre et al.l d2008l) . [Si/C] is expected to 
be ~ +0.7 indicating enrichment by Supernova e type II explo- 
sions. However, studying metal-rich systems ILevshakov et al.l 
(2009) recently find that these absorbers often have underabun- 
dant silicon with respect to carbon. The authors conclude that 
silicon is depleted into dust and, therefore, the absorber may by 
enriched by gas expelled from massive AGB-stars. 



7. Summary and conclusions 

We present a new, straightforward method to constrain the spec- 
trum of the intergalactic UV background by means of photoion- 
ization modeling of metal line systems. This new method re- 
quires simple, one-component systems with unblended, unsat- 
urated absorption features of various species. In particular, at 
least two different ionization stages of two different elements are 
needed. Then the column density ratio of two species of one ele- 
ment is used to construct a photoionization model of the system, 
while the second column density ratio is adopted to estimate the 
quality of the models for differ ent ionizing spectra. Base d on the 
UV background according to lHaardt & Madaul (1200 ll HM01) 
we parameterize three main characteristics of the spectrum: the 
Hen re-emission peak at 3Ryd, Ai, the intensity level in the 
range 3-4 Ryd, A2, which is related to the absorption due to the 
He n Lyman series, and the break at 4 Ryd, the ionization edge 
of Hen, A3. The three parameters are measured in terms of de- 
viation from the HM01 spectrum in dex. 

The method is tested with artificially constructed systems 
using Sim/Si iv to establish the photoionization models and 
Cm/Civ to evaluate the various SEDs. The main conclusions 
are: 

- The SED of the ionizing radiation cannot be reconstructed 
unambiguously. Instead, we find various spectra that repro- 
duce the observed column densities to a 1 cr confidence level. 

- The characteristics of the recovered SEDs are correlated. 
Evaluating the 1 cr confident models yields that Ai is cor- 
related with A2, i.e. the more pronounced the 3 Ryd Lya- 
peak the more depressed is the intensity between 3-4 Ryd. 
Therefore, we introduce (Ai - A2), which measures the 
break at 3 Ryd relative to the HM01 spectrum. Furthermore, 
(Ai - A2) is correlated with A3, which determines the hard- 



ness of the spectrum. This means, harder SEDs favor a more 
pronounced break at 3 Ryd. 

- Nevertheless, the main features of the ionizing spectrum can 
be recovered by estimating the average value of (Ai - A2) 
and inserting it into the (Ai - A2)-A3-relation. 

Moreover, we investigate the dependence of the estimated 
physical parameters of the metal line system, i.e. the ionization 
parameter log U, the metallicity [Si/H], and the relative abun- 
dance [Si/C], on the spectrum of the ionizing radiation. 

- The ionization parameter log U strongly depends on the 
3 Ryd Lya peak Ai, where higher A; leads to lower log U. 
This is due to adopting the Si m/Si iv ratio to constrain the 
photoionization model since the ionization potentials of Si m 
and Si iv are close to 3 Ryd (2.46 Ryd and 3.32 Ryd, respec- 
tively) making both column densities sensitive to the spec- 
trum in this energy range. 

- The metallicity correlates with the 4 Ryd break A3, confirm- 
ing the previous result that harder ionizing spectra lead to 
higher metallicities. 

- Considering all SEDs that reproduce the observed column 
densities within 1 cr confidence, the resulting values of the 
physical parameters are spread over Alogf/ ~ 0.5 dex, 
A[Si/H] ~ 1.5 dex, and A[Si/C] ~ 0.5 dex, respectively. 
Note, that in particular the uncertainty of the metallicity can 
be higher than one order of magnitude even for a well-fitting 
model if the spectrum of the ionizing radiation is uncon- 
firmed. 

However, the uncertainties decrease if the hardness of the 
SED can be constrained. For a fixed 4 Ryd break A3 the 
spread of log U and [Si/C] decrease by a factor of 2. In this 
case the spread of the metallicity decreases even more to 
A[Si/H] < 0.5 due to the strong dependence on the hardness 
of the ionizing spectrum. 

The method is applied to three appropriate metal line sys- 
tems at redshifts z ~ 1.75, 2.38, and 2.83. The recovered 
spectra differ from the HM01 background yielding the param- 
eters (Ai,A 2 ,A 3 ) = (+0.4, -0.5, +0.4), (+0.4,0.0,-0.3), and 
(+0.3, -0.6, -0.4), respectively. However, the recovered SEDs 
are subject to high uncertainties. 

(1) All systems favor an enhanced peak at 3 Ryd suggest- 
ing that the He 11 recombination re-emission is stronger than pre- 
dicted by the HM01 model. 

(2) At z < 2 the recovered SED is rather flat and hard. This 
is consistent with previ ous findings at similar redshifts by e.g. 
Agaf onova et al.l 00071) . However, the adopted column density 
ratios are Si n/Si m/Si iv and Al 11/ Al m, while for the z > 2 sys- 
tems we use Si m/Si iv and C in/C iv. Therefore, the hardening of 
the UV background at lower redshifts should be confirmed using 
UV data where the C in A911 line can be detected even at z < 2. 

(3) The recovered SEDs for the z > 2 systems tend to 
be softer than the HM01 background yielding a stronger break 
at 4 Ryd in combination with enhanced He 11 recombination at 
3 Ryd and intensity depression at 3-4 Ryd at least at z ~ 2.83. 
Ionization of the observed systems by a soft radiation field im- 
plies that the absorbers contain much He 11. Indeed, the expected 
He n/H 1 ratios are high, log rj ~ 2.4 and 2.8 at z ~ 2.38 and 2.83, 
respectively. This may suggest, that metal line systems reside 
in dense environments, where they are shielded from the gen- 
eral UV background and are exposed to a filtered radiation field 
(Mir alda-Escudell2005l ISchavdl2006l) . In particular, the system 
at z - 2.3799 towards HS 1700 +6416 is located wi thin a proto- 
cluster of galaxies reported by Steide l et al.l d2005l) . Generally, 
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metal line systems are thought to arise close to galaxies where 
the heavy elements are produced and then transported into the 
IGM by galactic winds. If this is true for the majority of metal 
line absorbers in QSO spectra, the usually adopted HM01 UV 
background might be too hard. Due to the strong dependence 
of the metallicity estimate on the hardness of the ionizing radi- 
ation, the derived metallicities might be overestimated by up to 
one order of magnit ude. However, the estimated metallicity of 
the diffuse IGM (e.g. lSchave e"taill2003t lAguirre et alj|2008l) is 
supposed to not be affected substantially. 

(4) Comparing the best-fit models adopting our recovered 
SEDs to models derived with theoretical UV background spec- 
tra from the literature we find a surprisingly well fit in case of 
the sa wtooth modulated UV background from Madau & Haardt 
(2009) assuming delayed helium reionization. Mad au & Haardtl 
d2009t) provide three different models boosting the He n/H i ra- 
tio by a factor of 10, 30, and 500 with respect to the standard 
case. The boosting factor of the best-fit model increases with 
increasing redshift. From direct observations of the He n absorp- 
tion it is known, that it shows a Lya forest structure and he- 
lium is therefore highly ionized at z S 2.7 dKriss et al.l [200 1 ; 
IZheng et alJ2004tlFechner et alJ2006atlShull et al.l20l6T) . Thus, 



the fraction of He n at these redshifts should be lower. However, 
the observed systems are probably exposed to a surprisingly soft 
radiation. They may reside in dense environments where remain- 
ing patches of He n might filter the UV background locally. 

The application of our new method to three observed sys- 
tems shows that it is possible in principle to infer constraints 
on the spectrum of the UV background. The recent finding of a 
rather soft ionizing radiation should be confirmed with a larger 
sample of systems in the future. However, systems exhibiting 
pronounced metal line features are probably residing close to 
galaxies. They may not be representative for the general inter- 
galactic medium since they live in special environments and may 
be affected by the radiation of local sources. Thus, they may be 
exposed to a softer radiation than the general UV background. 
As a preliminary result from the limit ed sample of systems in- 
vestigated so far, we suggest that the Ma dau & Haardtl |2009) 
spectrum with sawtooth modulation and delayed helium reion- 
ization is appropriate to stu dy metal line systems. In comp arison, 
the spectrum provided by iFaucher-Giguere et al.l d2009l) seems 
to be too flat and featureless although it might be a appropri- 
ate representation of the mean UV background in intergalactic 
space. 
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Fig. 13. Recovered SEDs in comparison to UV background 
spectra from the literature for the metal line systems at z — 
1.7529 (upper panel), z = 2.3799 (middle panel), and z = 
2.8266 (lower panel). Thick blue lines display our fiducial 
spectra described by (Ai,A2, A3) = (+0.4, -0.5, +0.4) at z = 
1.7529, (+0.4,0.0,-0.3) at z = 2.3799, and (+0.3,-0.6,-0.4) 
at z = 2.8266 , respe ctively, compared to the unmodified 
lHaardt & Madaul (UpOl) bac kground (thin black line; HM01), 
the sawtooth spectrum of iMadau & Haa rdt (200^1 in the 
lHaardt & Madaul ([1 996) uniform background limit (red dashed 
line; MH09/HM+S) and with delayed helium reionization 
boosting the Hen/Hi by factors of 10, 30, and 500 com- 
pared to the MH09/HM+S case (thin solid lines; MH09/DR10, 
MH09/DR30, MH0 9/DR500), respectively and th e UV back- 
ground spectrum of iFaucher-Giguere et al.l d2009l dotted line; 
FG09) at the appropriate redshifts. All spectra are normalized 
to log7 = -21.1 at 1 Ryd. 
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from the literature. Th e adopte d ionizing spectra are from 
lHaard t & Mad aul d2001l HM01) iFaucher-Giguere et all d2009l 
FG09), IMadau & Haardtl d2009l) assuming equilibrium con- 
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SEDs (see text). Shaded areas mark the observed ratio of 
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